DESY 93-139 
SCIPP 93/35 
November 1993 



Dijet Production at Large Rapidity Intervals 



Vittorio Del Duca 
Deutsches Elektronen-Synchrotron 
DESY, D-22607 Hamburg , GERMANY 

and 

Carl R. Schmidt Q 
Santa Cruz Institute for Particle Physics 
University of California, Santa Cruz, CA 95064, USA 



Abstract 

We examine dijet production at large rapidity intervals at Tevatron energies, by using 
the theory of Lipatov and collaborators which resums the leading powers of the rapidity 
interval. We analyze the growth of the Mueller-Navelet i^-factor in this context and find 
it to be negligible. However, we do find a considerable enhancement of jet production 
at large transverse momenta. In addition, we show that the correlation in transverse 
momentum and azimuthal angle of the tagging jets fades away as the rapidity interval 
is increased. 
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1 Introduction 



As the search for the top quark continues, the Tevatron Collider continues to produce 
large amounts of data on hadronic jets [l]. This data provides a unique opportunity to 
test our understanding of jet production at high energies. The calculation of jet pro- 
duction rates at hadron colliders is a challenging problem of perturbative QCD, because 
it involves many different scales, including: Aqcd, the hadron-hadron center-of-mass 
energy y/s, the parton-parton center-of-mass energy a/5, and the momentum transfer 
Q, which is of the order of the transverse momentum of the jets produced in the hard 
scattering. 

The conventional approach to these calculations is to work at fixed order in the 
coupling constant a s , assuming that V§, and Q are comparable in size, so that 
there are no large logarithms involving them. The effects of Aqcd are factorized into 
the parton structure functions, which are then evaluated at a scale of order Q using the 
usual DGLAP evolution. At present the first radiative corrections to the Born processes 
are available 0. These yield a more detailed description of the jet structure, reduce the 
dependence on the factorization scale, and are in very good agreement with the data on 
the one jet inclusive distribution at large transverse momenta [EJ. 

At the high energies of the Tevatron, however, there may be kinematic configurations 
where one cannot ignore the effects of the disparate energy scales. In the semihard region, 
defined as s » Q 2 >> Aq CD , the calculation of jet cross sections is characterized by 
the appearance of coefficients containing logarithms of large ratios of the kinematical 
invariants. If no restrictions are made on s = X\ X2 s, these logarithms will involve the 
small-x behavior of the structure functions, requiring a more sophisticated analysis than 
the usual DGLAP evolution (3J[. Combined with the experimental uncertainties in the 
structure functions at small-x, it appears very difficult to make precise predictions in 
this kinematic region. 

One way to overcome this problem is to try to disentangle the different ratios of 
kinematic invariants in the process. This can be achieved, for example, by requiring that 
the parton momentum fractions, X\ and X2, are large enough that no large ratios, other 
than the usual Q 2 /Aq CD , appear in the evolution of the parton distribution functions 0. 
The price to be paid is that the logarithms of the kind s/Q 2 will now appear in the parton 
subprocess. These logarithms are of the size of the rapidity interval in the scattering 
process. To realize this configuration experimentally, Mueller and Navelet proposed to 
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tag two jets at the extremes of the Lego plot in azimuthal angle and rapidity, at fixed xi, 
X2 and transverse momenta p±, and to watch the growth of the dijet inclusive cross section 
as the rapidity interval between the tagging jets grew with the center-of-mass energy. 
To deal with the large logarithms they used the Balitsky-Fadin-Kuraev-Lipatov theory 
(BFKL) ||, which systematically resums the leading powers in the rapidity interval by 
using a multigluon amplitude, with the gluons uniformly filling the rapidity interval 
between the tagging jets. The Mueller-Navelet i^-factor, defined as the ratio between 
the resummed and the Born dijet cross sections at large rapidity intervals and fixed 
x's, exhibits the power-like growth in the center-of-mass energy typical of the BFKL 
resummation. 

In this paper we study the Mueller-Navelet dijet cross section at the Tevatron energy 
of = 1.8 TeV. Since \fs is fixed, we instead let x\ and X2 vary with the rapidity 
interval. At the same time we retain the Mueller-Navelet requirement that xi, 22 are 
large enough that the parton distribution functions can be described by the DGLAP 
evolution. This is done by tagging on the two jets at the extremes of the rapidity 
interval y with transverse momenta larger than some cutoff p± m m- For reasonable values 
of y and p_i_ m in the momentum fractions x±,X2 will be sufficiently large. We can then 
study the effects of the minijets in the BFKL resummation as a function of the kinematic 
variables of the two tagging jets. 

In the exposition of this paper we follow the outline of Ref. @. Namely, in sect. 2 we 
consider the inclusive dijet production pp —> 2 jets + X at the Born-level, both exact 
H and in the large-rapidity limit. We then compute the leading logarithmic corrections 
at large rapidity as done in Ref. 0. In section 3 we present numerical results for the 
inclusive dijet production, with and without integrating over the jet transverse momenta. 
We discuss the effects of the BFKL resummation on the growth of the Mueller-Navelet 
K-factor, on the single jet pj_ distribution at large y, and on the correlations of the two 
jets in transverse momentum and azimuthal angle as the rapidity interval is increased. 
In section 4 we include some remarks on the range of validity and the limitations of this 
purely leading logarithmic calculation, and we present our conclusions. 

2 The Dijet Inclusive Cross section 

We are going to study the semi-inclusive process paPb — * 2 jets + X in the semihard 
regime defined by s » Q 2 , with Q 2 being a typical momentum scale in the event, 
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Q 2 ~ Pix P2L- The two tagged jets are chosen with a large rapidity interval y = y\ — y 2 ~ 
ln(s/pi^p2±)- Other relevant parameters in the event are the relative azimuthal angle 
and the rapidity boost y = (yi + 2/2) /2 of the two jets. 

In the semihard, large-y regime we can write the cross section: 

with = Q,Q,G labeling the structure function of the parton species and flavor = 
q, q, g inside hadron A(B). The parton subprocess cross section d&ij / dpf^dp^dcj) contains 
the sum over all additional particles (i.e. minijets) in the event. The factorization of the 
minijets into the subprocess cross section is possible, because at large y the initial parton 
momentum fractions x\ and x 2 are fixed in terms of the two tagged jet momenta, and 
are essentially independent of the particles filling the rapidity interval. We will arrive at 
this cross section in several steps, starting with the exact Born level cross section, taking 
it to the y » 1 limit, and finally filling in the rapidity interval with the minijets. 

i) Born Level Cross section . At the Born level the two partonic jets are produced 
back-to-back. The exact lowest order cross section can be put in the form (JJ) with the 
replacement 

M^^f '^-^) '(*-')■ W 

The parton momentum fractions and the subprocess invariants at this level are given by: 



x t = 


cosh(y/2) 

V s 


X2 = 


cosh(y/2) 

V s 


S = 


x 1 x 2 s = 2p 2 L (l + cosh(?/)) 


i = 


-p\{i + e-y) 


u = 





(3) 



where p± = pi± = p 2 ±. The subprocess invariants do not depend on the rapidity boost 
y. This is a general property, since y parametrizes the collective motion of the parton 
subprocess in the hadron reference frame. The lowest order parton cross sections are 
well known and can be found in, for instance, Ref. 
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ii) Large-?/ Born Cross section. We now investigate the lowest order cross section 
when the rapidity interval y is large. For y » 1 the lowest order amplitude is dominated 
by diagrams with gluon-exchange in the t-channel as in Fig. 1(a). In this limit the only 
subprocesses that contribute are gg — > gg and qg — > qg and qq — > qq. We obtain 

dt " 2p\ ' [V 
with Ca — N c — 3 the Casimir operator of the adjoint representation. Similarly, we find 



db qq Cf db qg Cp db gg 
di Ca di C\ dt 



(5) 



with Cp = (iV 2 — 1)/2N C = 4/3 the Casimir operator of the fundamental representation. 
Thus, it suffices to consider the subprocess gg — > gg and include the other subprocesses 
by means of the effective structure function 



f eff (x, /i 2 ) = G(x, /i 2 ) + ^ £[Q/(x, /i 2 ) + Q f (x, n 2 )] , (6) 

Ca f 



where the sum is over the quark flavors. The parton momentum fractions in the large-?/ 
limit are 

x = ^i e &+s/2) = —e Vl 

X2 = P^ e (-y + v/2) = P^ e -y, . (7) 

VS y/S 

Equation (0) is also valid in the large-jy limit when higher-order corrections are included, 
so that pij_ 7^ p2±- 

Minijet-corrected Cross section. As discussed in the introduction, going to higher 
orders in the coupling constant, i.e. to multiple parton emission, we encounter large log- 
arithmic contributions. In the semihard regime, the BFKL theory systematically 
resums the leading logarithmic terms \n(s/Q 2 ) by using a multigluon amplitude where 
the rapidity interval between the tagging jets is filled with gluons, strongly ordered 
in rapidity. This amplitude is shown in Fig. 1(b), where the thick line represents the 
resummation of the virtual radiative corrections, whose effect is to reggeize the gluons ex- 
changed in the t channel. The real gluons are inserted on these using the Lipatov effective 
three-gluon vertex [H . The BFKL multigluon amplitude is then put in a rapidity-ordered 
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phase space, the rapidities of the gluons are integrated out, and the dependence of the 
cross section on the gluon transverse momenta is reduced to the resolution of an integral 
equation. Its solution is then convoluted with a jet emission vertex on each side of the 
rapidity interval to give the minijet-corrected parton cross section for two jets at large-y: 



da. 



{19 



C A a s 



pI± 



f{y,Pi±,P2±) 



C A OLs 



pIx 



d 2 p 1± d 2 p 2 i_ 

In this equation f(y,Pi±,P2i_) is the Laplace transform in the rapidity interval y, 

f{y,Pi±,P2±) = J -^e" y fu(p 1± ,p2±), 
of the solution of the BFKL integral equation 



(8) 



fw(Pl±,P2±) 



(2vr) 2 



in(<j>— ir) 



-oo 



1XJ -V*H"(piJ 



-1/2-iv 



u — uiin, v) 



(9) 



(10) 



The eigenvalue of the integral equation u(n, v) is 

2C A a. 



uj{n, v) 



1 



7T 



Re^( 



\n\ + 1 



+ iv) 



(11) 



with ip the logarithmic derivative of the Gamma function. Substituting (^) and flTU| ) in 
(§), and doing the integral over u, the minijet-corrected parton cross section becomes 



da 



99 



cWs 



rOO / 

e W*-*) dve u ^ y cos lv In 



Pi± 

P2±, 



dpf ± dpl ± d(j) 4npl ± pl ± 
If we integrate over the azimuthal angle <ft i n fljjD, only the n = term survives. 



(12) 



iv) Minijet-corrected Cross section in the Saddle Point Approximation. At very large 
values of the rapidity interval y, the correlations between the two jets are washed out by 
the random walk in transverse momentum space of the gluons exchanged in the t chan- 
nel. This can be seen most easily by evaluating (12) in the saddle-point approximation. 
The contribution of fllT] ) to this equation is dominated by n = and is strongly peaked 
near v — 0. Thus we keep only the first term in the Fourier expansion in <f), and expand 
oj{y) = uj(0, v) about v — 



u (y) = A-Bv 2 + 



(13) 
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with 

A=^ln2, 5 = 1^C(3). (14) 

7T 7T 

Then we can evaluate ([12]) using the saddle-point approximation for the integral over is, 
to obtain 

da gg C\al e A y t ln 2 (p?±M^ X 



dp( ± dp 2 2L d(i) 8pi ± pj ± s/BTxy \ ABy 



eX P ~ ,o " J \ ■ ( 15 ) 



The exponential growth of ( p~5|) with the rapidity interval y is due to the production of 
the minijets. 



3 Numerical Results 

We now examine numerically the effects of the minijets at the Tevatron center-of-mass 
energy y/s = 1.8 TeV. We are mainly interested in understanding the behavior of the 
parton subprocess, which does not depend on y. Therefore, except where indicated, we 
work at fixed y and observe the cross sections as a function of the rapidity interval y. 
We chose y = so that neither x can become too small. For consistency of notation 
we will refer to the leading jet in rapidity as jet 1 and the trailing jet as jet 2 (i.e., 
Vi — +M/2, yi = — M/2). Of course, everything is symmetric under the exchange 



of the two jets. We have used the leading order CTEQ structure functions ]ItJ with 
the renormalization and factorization scale set to the geometric mean of the transverse 
momenta of the tagging jets, fi 2 = pi±p2j_- We shall address some of the difficulties 
involved in the choice of scale in the next section. For p_i_ m m > 10 GeV and y = the 
parton density functions are always evaluated at x > 10~ 2 , so we are justified in using 
the DGLAP evolution in this region of phase space. 

We begin by looking for the exponential growth of equation (|To| ) in the cross section as 
originally suggested by Mueller and Navelet. To do this we integrate over the azimuthal 
angle <fi and over both transverse momenta above a cutoff of p_i_ m i n = 20 GeV. In Fig. 2(a) 
we present this cross section in the first three approximations (i-iii) of section 2. From 
the plot we see that the large-y Born cross section is a good approximation to the exact 
Born level cross section for large y ^ 4. However, the minijet-corrected cross section 
does not exhibit any great enhancement at large rapidity. This is more easily seen in a 
plot of the K-factor, defined here as the ratio of the minijet-corrected cross section to 
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the large-y Born cross section, 



rfq (minijet) / do-(large-y) 
dydy I dydy 

The i^-factor is defined so that K — ► 1 as y — > 0. In Fig. 2(b) we see that the .ff -factor 
increases until y ~ 6, but then quickly goes to zero Q 

This effect can be understood if we remember that the rapidity dependence enters 
not only in the BFKL kernel f(y,Pi±,P2±), but also in the parton structure functions 
fi(xi, fj J 2 )fj(x2, yit 2 ) where the momentum fractions are given by (|7|). The allowed phase 
space in pi_i_,p2_i_ is substantially decreased at large y by the restriction that the momen- 
tum fractions must be less than 1. The decrease in phase space has a greater effect on the 
minijet cross section, with the result that the i^-factor is cut off at large y. It requires 
a much larger range in scales from p± m { n to */s in order to approach the exponential 
growth of the Mueller-Navelet .fT-factor. This is exhibited in Fig. 3 where we show the 
.fT-factor at the Tevatron energy, at a Large Hadron Collider energy of \fs = 15 TeV, 
and at y/s = 10 5 TeV. The full exponential growth is achieved in the limiting case of 
an infinitely large value of \/s, where the cutoff in the phase space of the minijet cross 
section never occurs. 

Thus, at the Tevatron energy we must look elsewhere for effects of the minijets. In 
Fig. 4(a) we show the p± distributions of jet 1 for a rapidity interval of y = 4. We plot 
the minijet cross section with two different cutoffs for the second jet, P2_i_min = 10 GeV 
and 20 GeV, while the Born level cross section always has p2± = pi±- These plots exhibit 
two effects of the minijets. We see that the overall scale of the p± distribution depends 
strongly on the minimum p± of the second jet, and that the slope of the distribution 
is flatter with a substantial increase at large pj_. In Fig. 4(b) at y = 6 we see an even 
greater dependence on P2±min- These effects can be partially understood by the fact that 
the p± of jet 1 can be balanced by the smaller p± of jet 2 and the minijets produced 
in the rapidity interval. The lower x values required for this type of event increase its 
likelihood relative to the back-to-back dijet event, which is all that can occur at the Born 
level. It is even possible for jet 1 to attain transverse momenta that are kinematically 
impossible at the Born level. 

These arguments suggest that the minijets occuring in the rapidity interval between 
the tagged jets will cause the tagged jets to become uncorrelated. This can be seen easily 

2 The fact that K goes below 1 for small y is presumably a computational artifact, arising from the 
difficulty in doing the numerical integration over a very sharply peaked function at small y. 
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in the minijet formulae by looking at the BFKL kernel as it is varied from y = to very 
large y. For very small rapidities we approach the Born cross section with 



f(y,Pi±,P2±) 



(17) 



The two jets are produced back-to-back in p± and 0. However, as the rapidity interval 
becomes large, we have 



and the tagging jets become completely uncorrected. 

The disappearence of correlations as y increases can be seen dramatically in Fig. 5(a) 
where we plot the transverse momentum distribution of jet 1 at a fixed value of p2± = 50 
GeV. For a rapidity interval of y = 2 the cross section is strongly peaked near px± = p 2 ±. 
As the rapidity is increased there is a diffusion of the jet 1 momentum away from the 
jet 2 momentum until the peak is practically gone for y = 5. In practice jet 2 will be 
integrated over some range of transverse momenta, so in Fig. 5(b) we show the same plot 
with p2± integrated from 50 GeV to 55 GeV and \y\ < 0.5. To retain the normalization 
we have divided this cross section by 5 GeV. 

Similarly, there is also a reduction of the correlation in the azimuthal angle as the 
rapidity interval increases. This can be seen in Fig. 6, where we show the distribution 
with both jets integrated from p_i_ m m = 20 GeV. The distribution is normalized to the 
uncorrelated cross section da/dydy, so that the area under each curve is equal to 1. As 
expected the correlation in decreases as we vary from y = 5 to y — 7. The decorrelation 
in 0, however, is slower than the decorrelation in p±, because the eigenvalue (|TT|) of the 
BFKL integral equation is more strongly peaked in v than in n. For example, at y — 5 
the tagging jets are not correlated any more in p±, while they still show a considerable 
correlation in 0. 

4 Discussion and Conclusions 

The BFKL analysis that we have been using is a leading logarithmic approximation. 
With this in mind we offer some caveats to our results and discuss which effects should 
survive in an exact calculation. First, we must state that any of our plots at very 
small y are not expected to be very accurate. For y ^ 2 there is even a reasonable 




(18) 
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discrepancy between the large-?/ Born Cross section (ii) and the exact Born level cross 
section (i). However, we expect that the trends as y is increased should be apparent 
even at reasonably small values of the rapidity interval. In particular the decorrelations 
in p± and <fi should definitely increase with y. 

There are also several ambiguities in our calculation, arising from the fact that the 
BFKL analysis assumes little variation in the p± of the minijets. For instance, the 
rapidity variable used in the standard BFKL analysis is Y = ln(s/Q 2 ) where Q is some 
typical scale of the minijets. In our calculation we have chosen Q 2 = pi±p 2 ± so that 
Y = y, the experimental rapidity. This should not make a significant difference at large 
y, but it does emphasize the fact that the approximation becomes less reliable when 
the transverse momenta of the tagging jets are not of similar size, and there arise large 
logarithms of order ln(pi±/p 2 ±). In addition, the dependence on jet definition, cone size, 
and other variables are subleading in this analysis at large y. 

For related reasons the proper renormalization/factorization scale fi at which the 
coupling constant is evaluated is not well-determined in our leading logarithmic analysis. 
As required for the BFKL solution, we have evaluated the coupling for all of the minijets 
at a single scale of order Q. Possible choices for /i 2 are pi_i_p2_i_ (as we have used), p 2 ± , pl ± , 
or max(p 2 _|_, p 2 L ) . At the level of our approximation, all of these scales are equivalent, 
but in practice the choice of scale can make a reasonable difference in the slope of the 
Pj_ spectrum []. However, our main conclusions about the decorrelation in transverse 
momentum and azimuthal angle at large y, as well as the increase in the cross section 
at large p± and y will not change. 

Finally, it is interesting to imagine a comparison of our results with a fixed 0(a^) 
calculation. The fixed order calculation includes only the effects of up to three parton 
jets, while the BFKL resummation includes the leading contributions of amplitudes 
containing an arbitrary number of parton jets. From naive estimates one might expect 
that y ^ 6 is not large enough to warrant the use of the full BFKL analysis and that 
a next-to-leading order calculation is quite sufficient. However, as we have seen, the 
kinematic phase space is greatly enhanced by the sharing of transverse momentum among 
the additional minijets. This suggests that it may be necessary to go beyond 0(a^) at 
large rapidity intervals. Moreover, the BFKL approximation clearly predicts the main 
features of the multiple jet emission at large y and readily suggests the experiments to 
look for them. It would be exciting to compare the predictions of the BFKL resummation 

3 This ambiguity will arise in any non-exact calculation involving multiple jets with disparate energies. 
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given here, as well as the next-to- leading order calculations, against experiment. 
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Figure captions 

Fig.l: Two jet production amplitude in the large-y limit at (a) the Born level and (b) 
with minijet corrections. 

Fig. 2: (a) Inclusive dijet production at the Tevatron, as a function of the rapidity interval 
y. The dashed and dot-dashed lines are respectively the exact and large-?/ Born 
cross sections, and the solid line is the minijet-corrected cross section, (b) The 
-fT-factor, i.e. the ratio of the minijet-corrected cross section to the large-?/ Born 
cross section, as a function of the rapidity interval y. The kinematic parameters 
for both figures are described in the text. 

Fig. 3: i^-factor as a function of the rapidity interval y at different center-of-mass energies. 
The cutoff on integration over both transverse momenta is at p_i_ m in = 20 GeV. From 
bottom to top, the solid lines represent the i^-factors at Tevatron energies, at LHC 
energies (y/s= 15 TeV), and at y/s = 10 5 TeV. 

Fig.4: p± distribution of jet 1 at (a) \y\ = 4 and (b) \y\ = 6. The dashed and dotdashed 
lines are respectively the p± distributions for the exact and the large-y Born cross 
section, for which p 1± = p 2 _i_. The solid lines are the p± distributions of jet 1 for the 
minijet-corrected cross section, with two different cutoffs for jet 2, p 2 ±min = 10 GeV 
and 20 GeV. Notice that in (b) the dashed and dotdashed lines completely overlap. 

Fig.5: distribution of jet 1 with the jet 2 transverse momentum (a) fixed at 50 GeV, 
and (b) integrated from 50 GeV to 55 GeV. From top to bottom, the solid lines 
are the p± distributions for the minijet-corrected cross section at \y\ = 2, 3, 4, 5 
and 6. In (b) the rapidity boost y is integrated over \y\ < 0.5. 

Fig.6: distribution normalized to the uncorrelated cross section dajdydy. From top 
to bottom, relative to the peak, the solid lines are the distributions for the 
minijet-corrected cross section at \y\ = 5, 6 and 7. 
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